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We demonstrate that images of flux vortices in a superconductor taken with a transmission elec¬ 
tron microscope can be used to measure the penetration depth and coherence length in all directions 
at the same temperature and magnetic field. This is particularly useful for MgB 2 , where these quan¬ 
tities vary with the applied magnetic field and values are difficult to obtain at low field or in the c 
direction. We obtained images of flux vortices from a MgB 2 single crystal cut in the ac plane by 
focussed ion beam milling and tilted to 45° with respect to the electron beam about the crystallo¬ 
graphic a axis. A new method was developed to simulate these images which accounted for vortices 
with a non-zero core in a thin, anisotropic superconductor and a simplex algorithm was used to make 
a quantitative comparison between the images and simulations to measure the penetration depths 
and coherence lengths. This gave penetration depths Aab = 100 ± 35 nm and Ac = 120 ± 15 nm 
at 10.8 K in a field of 4.8 mT. The large error in Aab is a consequence of tilting the sample about 
a and had it been tilted about c, the errors on Aab and Ac would be reversed. Thus, obtaining 
the most precise values requires taking images of the flux lattice with the sample tilted in more 
than one direction. In a previous paper [J. C. Loudon et al., Phys. Rev. B 87, 144515, 2013], we 
obtained a more precise value for Aab using a sample cut in the ab plane. Using this value gives 
Aab = 107 ± 8 nm, Ac = 120 ± 15 nm, ^ab = 39 ± 11 nm and ^c = 35 ± 10 nm which agree well 
with measurements made using other techniques. The experiment required two days to conduct and 
does not require large-scale facilities. It was performed on a very small sample: 30 x 15 pm and 
200 nm thick so this method could prove useful for superconductors where only small single crystals 
are available, as is the case for some iron-based superconductors. 

PACS numbers: 74.25.Uv, 68.37.Lp, 74.25.Ha, 74.70.Ad 
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I. INTRODUCTION 

Superconductors have zero electrical resistance and ex¬ 
pel magnetic flux from their interiors (the Meissner ef¬ 
fect). However, if a sufficiently high magnetic held is 
applied, hux penetrates by Rowing along channels called 
Rux vortices. Each vortex carries one quantum of mag¬ 
netic flux, $0 = /i/2e where h is Planck’s constant and e 
the electron charge. They consist of a core with a radius 
given by the coherence length, where the number of 
carriers (electrons or holes) contributing to superconduc¬ 
tivity is suppressed. Electrical supercurrents circulate 
around the centre, diminishing over a radius given by the 
penetration depth, A. In a conventional superconductor, 
the coherence length is related to the energy required to 
excite a carrier out of the superconducting state. A, and 
the velocity of the carriers at the Fermi energy, vp, via 
^ = hvp/T^^ and the penetration depth is related to the 
number density of carriers involved in superconductivity, 
ns, and their effective mass, to*, via A = m*/pL^nse^ 
(/To is the permeability of free space). 

In a type-I superconductor, the core exceeds the size 
over which the supercurrents persist and vortices attract 


one another as the area of normal (non-superconducting) 
material is minimised if the cores overlap. In a type-II 
superconductor, the supercurrents persist over a larger 
radius than the core and the Lorentz force causes vortices 
to repel one another so they form a hexagonal array in 
an isotropic superconductor. Introducing the Ginzburg- 
Landau parameter, k = A/^: a type-I superconductor 
has K < and type-II has k > I/-\/2- 

An anisotropic superconductor has different properties 
along different crystal axes, a, b and c. Most are uniax¬ 
ial so that a and b are equivalent. The anisotropy in 
the penetration depth is 7 a = Ac/Aab and in the coher¬ 
ence length it is yj = U&/fc- In a 1-band superconduc¬ 
tor, where there is one source of carriers contributing to 
superconductivity, the penetration depth and coherence 
length are independent of the applied magnetic field and 
their anisotropies are equal. One method to investigate 
the penetration depths and coherence lengths in both the 
a and c directions is to induce flux vortices with their axes 
normal to the ac plane. The vortex then has an elliptical 
core surrounded by circulating currents as illustrated in 
Fig. [Ha). 

MgB 2 is a rare 2-band superconductor— discovered in 
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FIG. 1. (a) Flux vortex with its axis normal to the ac plane. 
The core (shaded) is elliptical with dimensions ^ab and ^c- 
Supercurrents J follow ellipses with the same aspect ratio, di¬ 
minishing on lengthscales Ac and Aab- Streamlines of J are 
also contours of magnetic flux density B. (b) Experimental 
arrangement for imaging flux vortices. The electron beam 
is deflected by the component of B normal to the beam so 
vortices appear as black-white features in an out-of-focus im¬ 
age. (c)-(f) Simulated images with defocus A/ = 7.15 mm 
for a flux vortex in a 180 nm thick specimen in the orien¬ 
tation shown in (a) but tilted 45° about a. Contrast values 
(see text) are shown below each image, (c) Aab = 100 nm, 
Ac = 120 nm, ~ 34 nm, (d) Aab = 100 nm, Ac = 120 nm, 
= 1 nm, (e) Aab = 200 nm, Ac = 120 nm, = 34 nm, (f) 
Aa6 = 100 nm, Ac = 200 nm, = 34 nm. 


2001 with a transition temperature = 39 K. It is uniax¬ 
ial with a hexagonal crystal structure^ (space group 191: 
P6/ mmm) composed of alternating layers of magnesium 
and boron with lattice parameters a = b = 3.086 A and 
c = 3.542 A. Two bands contribute to superconductiv¬ 
ity: the (T-band associated with bonding from the boron 
Pxy orbitals and the 7r-band associated with boron pz or¬ 
bitals^. a carriers are confined to the ab planes whereas 
the TT carriers are delocalised almost isotropically. At low 
magnetic fields, both cr and tt bands contribute to super¬ 
conductivity but as the field is increased, the tt contribu¬ 
tion diminishes so that above 0.8 T (at 2 K), only the a 
band contributes^. This has the effect that the penetra¬ 
tion depth and coherence length vary with field^. 

In a previous paper—, we showed that the penetration 
depth of MgB 2 , Aab, could be obtained in the low-field 
limit by making a quantitative comparison between im¬ 
ages of flux vortices acquired using transmission elec¬ 
tron microscopy and simulations. Here, we extend this 
method and show that the penetration depths Aab and 
Ac and coherence lengths ^ab and ^c can be measured in 
a low field of 4.8 mT from a very small sample. 

Focussed ion beam milling was used to cut a MgB 2 
sample in the ac plane of size 30 x 15 pm, thinned to 
200 nm so that it was electron transparent (see sec¬ 
tion nni). Flux vortices penetrate normal to the thin 
surfaces and the sample was tilted about its a axis at 


a = 45 ± 5° to give a component of the B-field normal 
to the electron beam (Fig. HKb)). The beam is deflected 
by the Lorentz force and flux vortices appear as black- 
white features in an out-of-focus imaged. Such images 
are sensitive to the H-field throughout the thickness of 
the specimen, not just the stray field, and they can be 
acquired in real time^. 

The effect on such images of changing the coherence 
length and penetration depth is shown in Fig. [Ilc)-(f). 
These images were simulated by extending Beleggia’s 
method^ to model vortices with a non-zero core in a 
thin, anisotropic superconductor (see section |TT1) • We use 
Klemm and Clem’s Ginzburg-Landau model for the vor¬ 
tex core^iS. In this model, the core has the same symme¬ 
try as the circulating currents so that ^ab/^c = ^c/^ab 
and = Cab^^c^^ but any model for the magnetic struc¬ 
ture of a flux vortex could be used with equal facility. 

Fig. He) shows a simulated image of a vortex with 
Aab = 100 nm, A^ = 120 nm and = 34 nm. (d) 
shows that decreasing to 1 nm sharpens the image, in¬ 
creasing the contrast (the difference in the maximum and 
minimum intensities divided by their sum) from 10.5% to 
14.8%. In (e), Aab is doubled which stretches the image 
in c and reduces the contrast from 10.5% to 8.8%. (f) 
shows the images are most sensitive to A^ so that when 
Ac is increased to 200 nm, the image is stretched in a and 
its contrast falls to 4.4%. This sensitivity of the images 
to changes in these parameters should allow the simulta¬ 
neous measurement of the penetration depths and coher¬ 
ence lengths in all directions. In this paper we assess the 
accuracy of this new technique. 


II. SIMULATION OF FLUX VORTEX IMAGES 

In this section, we present a model to calculate the 
magnetic fields generated by a flux vortex and from this 
simulate transmission electron micrographs. The model 
accounts not only for the B-field inside the supercon¬ 
ductor but also for the spreading of the field lines near 
the superconductor surface and the field outside. It ex¬ 
tends Beleggia’s method^dS to treat the case of a vortex 
in a thin, anisotropic superconducting slab and makes 
use of the work of Klemm and Clent^i^iii to account for 
a non-zero vortex core although it has the convenient 
feature that any model for the vortex core can be used 
with equal facility. Like all magnetic objects, flux vor¬ 
tices change only the phase and not the intensity of the 
electron beam and once the fields have been calculated, 
the phase shift can be found using the Aharanov-Bohm 
formula. Once the phase shift is known, any image can 
be simulated. Thus we first evaluate the magnetic vector 
potential, then use this to find the phase shift and from 
this simulate out-of-focus images of flux vortices. 
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FIG. 2. The relationship between the coordinates X, Y, Z, re¬ 
ferring to the specimen and the microscope coordinates x,y, z. 
The y and Y axes are normal to the other two axes and point 
in the direction given by the right-hand rule. The axis of the 
vortex is parallel to Z and the B-field enters the specimen at 
the bottom and exits at the top. The electron beam is parallel 
to z. 


If the order parameter of the superconducting state is 
written in terms of its amplitude / and phase d as d' = 
/e*®, the vector potential inside the superconductor is 
related to it by the 2nd Ginzburg-Landau equatio n^^d^ : 

Af -f LV X (V X A) = p-fX9 (2) 

27T 

where L is penetration depth tensor. When X, Y, Z are 
principal axes of the superconductor it has components: 

/A2 0 0 \ 

A = 0 A2 0 (3) 

V 0 0 A^J 

Following Clem^, we look for a solution of the form 


A. Coordinate Systems 


A — Abuik T ^surface 


( 4 ) 


In order to visualise flux vortices using transmission 
electron microscopy, the specimen must be tilted by an 
angle a to provide a component of the B-field normal to 
the electron beam so that the electrons are deflected by 
the Lorentz force and show contrast in an out-of-focus im¬ 
age. Thus we follow Beleggia’s method^i^ and introduce 
two coordinate systems: A, Y, Z referring to the speci¬ 
men with X and Y in the specimen plane and Z normal 
to its surface and x, y, z referring to the microscope with 
z parallel to the electron beam as illustrated in Fig. [2] 
The specimen surfaces are aX Z = zLd so its thickness is 
t = 2d. We first evaluate the magnetic vector potential 
in terms of the specimen coordinates and use this to find 
the phase shift in the xy plane which is equivalent to the 
plane in which images are recorded. 


Abuik is the solution for a single vortex in an infinitely 
thick specimen and consequently has no ^-dependence. 
Asurface is the general solution with the correct boundary 
conditions but with the right-hand side of the Ginzburg- 
Landau equation set to zero. 

Cleroi^ solved for the bulk term in the isotropic case 
by using an order parameter for a single vortex of the 
form 


= /e*® 




( 5 ) 


where p is the radius from the axis of the vortex, (j) 
is the azimuthal angle and ^ is the coherence length. 
Klemm and Glerr^di later extended this solution to the 
anisotropic case so that the order parameter becomes: 


B. Magnetic Vector Potential 

Here we evaluate the magnetic vector potential, A, 
from a flux vortex passing through a thin supercon¬ 
ducting slab with its axis directed along Z, normal to 
its surfaces. We use Beleggia’s Fourier-space methodic 
throughout in which all the functions are Fourier trans¬ 
formed in X and y but not z (or X and Y but not Z). 
This allows the Fourier transform of the vector poten¬ 
tial and phase to be expressed by analytical but very 
lengthy expressions. These were evaluated symbolically 
using Matlab and only the final inverse transforms were 
performed numerically. For both coordinate systems, we 
use the transform convention that if g(x,y,z) is a func¬ 
tion in real-space, its Fourier transform g(fcx, ky, z) is: 


g{kx,ky,z) 


J J 5(x, y, dxdy (1) 




I AVA2+r2/A| 

A2/A2 +y2/A2 +^2/A2 


o-i arg(X/A +iY/K ) 


( 6 ) 

where A = (A^ A^. This gives the magnetic flux 

density as 


P Ko{R) ^ 

2^A,A, (^/A)Ai(e/A) ^ ^ 

where R = (A/Ay)^ -|- (F/A^)^ -I- (^/A)^ and Aq 

and Ki are zero and first order modified Bessel functions. 
Fourier transforming the flux density gives 


$oAi(Qe/A) 9 

QKi{i/X) 


where Q = y^I -|- A:2 A2 -|- fc^A^. Since B = V x 
in Fourier space we have B = (ik^Ay — ik.^Ax)'Z 


( 8 ) 

A, 


— oo —oo 
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X A where kj^ = {k^ ,ky,Q). Imposing the additional 
requirement that A obey the London gauge V.A = 0 so 
that ik^Ax + iky Ay = 0, the vector potential in Fourier 
space is 


klQK,{^/A) 



(9) 


We now find the surface term, Agurface- This is the 
general solution to the Ginzburg-Landau equation but 
where the right-hand side is set equal to zero. We use 
the approximation introduced by Clem^ that the surface 
term should be only weakly influenced by the vortex core 
and so set /^ = 1 in the search for a solution. The 
validity of this simplification was confirmed by Brandli^ 
who modelled the complete vortex and found that the 
core only expands by a few percent as it approaches the 

I 


surface of the superconductor. The surface term is thus 
the solution to 


A-t-TV X (V X A) = 0 (10) 

It should be noted that setting the right hand side of 
the Ginzburg-Landau equation (eqn. [2]) to zero fixes the 
gauge of the vector potential and for an anisotropic su¬ 
perconductor, this is not the London gauge. Thus we 
cannot say V.A = 0 nor make the convenient replace¬ 
ment V X (V X A) = — V^A. Instead we must deal with 
the awkward cross-terms arising from the double curl. 

Taking the Fourier transform of Eqn. 1101 gives: 


A -I- L{i kj^ -I- Z d^) X ((i kj_ -I- Z 9^) X A) = 0 (11) 

We now postulate a solution of the form 
A{k^ ,ky,Z) = a{k^ , ky)e^^ and the resulting equation 
can be written in matrix form as: 


/ 1 + klAl-AlP^ 

-k^kyAl 


-k^kyAl 
l + fc 2 A 2 _a2/32 

ikyA^P 


ik^A^ P 
ikyA^P 

1 +A|(fci +kl) 



( 12 ) 


The above equation gives non-zero solutions for the 
vector potential if the matrix cannot be inverted. To 
achieve this, values of P must be found to make the de¬ 
terminant zero. At this point we introduce the symme¬ 
try of the problem otherwise the answers become very 
lengthy. For this experiment, the specimen was tilted 
about the a axis and thus, A || c, T || a and Z \\ b so 
that = Ac and A^ = Aaf, and A^ = Aab- There are 
then four possible values of /3: 

/3i.3 = ±Qa/Aab (13) 

and 

P2,4 = ±Q/Ac (14) 


where Qa = \J^ + {k\+ kl)Al^ and Q = 

+ k'j^Aab + kyAc. The complete vector potential in¬ 
side the superconductor is then: 

4 

Ajnside — Abulk T ^ (1^) 

n —1 

where Ci - C 4 need to be determined by the boundary 
conditions. 

This leaves the vector potential outside the super¬ 
conductor to be determined. Maxwell’s third equation 
gives V X B = 0 as there are no electrical currents out¬ 
side the superconductor so the vector potential obeys 
V X (V X A) =0. This time, the London gauge, V.A = 0, 
may safely be invoked to give 


with corresponding eigenvectors 


0 


— I ^iQa/i^ky Adb') 


and 


(15) 


V^A. 


outside 


= 0 


or, in Fourier space: 


, 2 T 


= 0 


(18) 


(19) 


The solution to this is 


±^AS + klAl,)/{Al,k,Q) 
a2,4 = I ±.1 Achy IQ 

1 


(16) 


Atop — 


C 5 

C 6 I e 

i{c^k^ + c^ky)lky 


-k,Z 


( 20 ) 
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Abottom = j cs j (21) 

\ -i{c7k„ + C 8 fcy)/fc_L ) 

where the Z-components are determined by the Lon¬ 
don gauge. 

We can now simplify the equations as symmetry re¬ 
quires that Ay,y{—Z) = and A^{—Z) = 

—A^{Z). This gives ci = —C 3 , C 2 = —C 4 , C 5 = c-j and 
ce = cs- 

Summarising so far, the vector potential inside the su¬ 
perconductor is 


must be continuous across the boundaries at Z = ±d as 
there are no currents confined to the surface of the super¬ 
conductor. There is also the requirement that the normal 
component of the flux density B± be continuous but this 
arises from Maxwell’s 3rd equation, V.B = 0, and by 
using a vector potential, it is automatically satished. 

Condition (1) that the in-plane vector potential is con¬ 
tinuous at Z = zLd gives two equations (one for each 
component): 

Abuik.x + 2ciai ^ cosh(/3id) -I- 20203 ^ cosh(/ 32 d) 

= 056 "'=^"* (24) 


( 0 ^ 3 ,. cosh(/?iZ) \ 

Oi y cosh(/3i.^) I 

Oi ^ sinh(/?iZ) J 

/ cosh(/32^) \ 

-I- 2 c2 I Y cosh(/32^) (22) 

y 03 ^ sinh(/? 2 ^) j 

and above and below the superconductor, the vector 
potential is 

/ C5 \ 

Aabove, below =1 Og I eA ^ (23) 

y ±*( 05 ^ 3 ^ +C(,ky)/k_L ) 

To fix the values of Ci, C 2 , C 5 and ce, we invoke the 
following boundary conditions: (1) In order to calculate 
the phase shift from the vector potential, the X and Y 
components of the vector potential must change contin¬ 
uously across the boundary between the superconductor 
and vacuum at Z = id. (2) The in-plane flux density i?|| 

I 


Auik.v + y cosh(/3id) i 20203 cosh(/ 32 (i) 

= (25) 

Calculating the flux density via B = V x A or, in 
Fourier space, B = (ikj^ izcl^) x A and matching its 
in-plane components at the interface gives two more: 

{-(k^/ki_){c^k^ i ceky) + cekj_)e~''^‘^ = 

2ciDi sinh(/3i(i) i 2 C 2 D 2 sinh(/ 32 d) (26) 


(-csfcj. i {k^/k±){c 5 k^ i ceky))e = 

2 C 1 D 3 sinh(/3i(i) i 2c2i44 sinh(/32d) (27) 

where Di = — Pia^ y, D 2 = ik^a^ ^ — P 2 a 2 ,Y^ 

D 3 = /SiOj 3 ,. — ik^a^ z and D 4 = P 2 a 2 ,x ~ 

Writing these four equations in matrix form gives: 


/ 


2 II 1 sinh(/3id) 2 II 2 sinh(/ 32 d) 


2113 sinh(/3id) 2 II 4 sinh(/32d) 

2ai,x cosh(/3id) 202 , 3 ; cosh(/l 2 d) 
y 2oi,y cosh(/3id) 2 o 2 ,j/ cosh(/ 32 d) 


„ — k I d 

C- 


'^-k±\e 

_g—fc I'd 

0 


— k±d 


^-k± 


^ — k\ d 


^x ^Y „ — k I d 

k± ® 

0 

„ —A; I d 


/ Cl \ 

( ' \ 

C. _ 

0 

C 5 

_ A 

■^hulk,X 

V C6 / 

\ “Aulk.V / 


(28) 


The coefficients ci, C2, C5, ce can then be found by 
inverting the matrix and then, after substituting the an¬ 
swers into Eqns. [22] and |23l the vector potential is fully 
determined. It should be noted that although we use 
Klemm and Clem’s solution for Abuik; Eqn.|5S|shows that 
our method has the convenient feature that any model for 
Abuik could be used with equal facility. 


C. Phase Shift 

The magnetic contribution to the phase shift suffered 
by the electron beam after passing through a specimen is 
related to the vector potential via the Aharanov-Bohm 
expression: 


A(x, y, 2 :).dl 


Hx,y) 


27re 

h 


—00 


(29) 
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1 = (X — w tan a)! + FJ + wK. 


(33) 


By differentiating the above equation, an increment in 
its trajectory dl can be related to an increment in w via: 


dl = (— tan a I + K)di 


(34) 


FIG. 3. The relationship between the coordinates X, Y, Z, 
referring to the specimen (shown by the grey rectangle) and 
the microscope coordinates x,y,z. 


The phase shift written in terms of the specimen coor¬ 
dinates is thus: 


4,{X,Y) = -‘^x 


— tana 


/ OO / — LCLil LX. \ 

A{X — wtana,Y,w). ( ^ j 

In the previous section the vector potential was cal- 
and those in X, Y and Z of the specimen coordinate culated in Fourier space. This is related to the vector 
system are denoted I, J and K, the two sets of unit potential in real-space via the inverse Fourier transform: 


where dl is an increment along the trajectory of the 
electrons shown in Fig. [3] If unit vectors in x, y and z of 
the microscope coordinate system are denoted i, j and k 


vectors are related via: 


A{X,Y,Z) = 


OO OO 


I = i cos a -I- k sin a 

J = -j 

K = i sin a — k cos a 


(30) 

(31) 

(32) 


^ J J A{k,,k,,Z)e^^^-^dk,dk, (36) 

— OO —OO 

where R = (2f, Y^Z). Thus it follows that 


If the electron passes through the specimen at position 
(X, y, 0) and if for the preceding part of its journey, we 
label the height above the specimen in the Z direction w 
(see Fig. El, its position at any point in its trajectory is 
given by 


A(X — wtana, Y, Z) = 

ff dfc^(37) 


So the phase shift is 


27re 


^{X,Y) = -— I ( — // A(fc,,fc,,«;)e*'^-^e-*'=-’"‘““dfc,dfc, 


— tana 
0 
1 


din 


(38) 


It can be seen from Fig. [3] that if an electron passes through a point (A, Y, 0) on the specimen, it passed through a 
point (x,2/,0) in the xy plane where x = X cos a and y = —Y. The phase shift in the xy-plane (which is equivalent 
to the plane in which the image is taken) is now found by making this substitution. 


o-—p poo / ^ r r \ / tan a \ 

cl){x,y) = - —J I q 


(39) 


Now let kx = k^/ cos a and ky = —ky-. 
27re f°° / 1 


(j)[x, y) = 


/ 


h 7-00 \47r2 


Aikx cos a, —ky, e 


sin a ^ 


— tana 
0 
1 


dw (40) 


Changing the order of integration gives: 
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(t>{x,y) = 


2 -Ke cos a 1 
h 47r^ 


^'ik_L .r 


A{kx cos a, —ky, w)e 


—ikxW sin a 


die 


\ J 

1 — tan a A 

j dkx dky. 

0 

1 / 


(41) 


or representing an inverse Fourier transform as IFT and using the flux quantum 4>o we have: 



(j){x,y) = IFT 


-400 -200 0 200 

jc/nm 


. — Sin a 
° ' cos Of 


-400 -200 0 200 

x/nm 


A{kx cos a, —ky, w)e 


— iks 


w sin a 


dw 



(42) 


I- 

be taken. Taking an out-of-focus image is equivalent to 
propagating the wavefunction through free space by a 
distance A/, known as the defocus. This is done via the 
Fresnel-Kirchoff integral^ so that the defocussed wave- 
function tpAfix, y) is related to the in-focus wavefunction 
via: 


ti)Af{x,y) = 
1 


AA/ 


JjMx' dx'd2/'(43) 


FIG. 4. Contours of the phase-shift spaced by 0.2 rad illus¬ 
trating the projected B-field from a flux vortex with Kab = 
100 nm and Ac = 120 nm and (a) = 1 nm and (b) 

= 36 nm. The red lines indicate the surfaces of the speci¬ 
men which is 128 nm thick. The a-axis is parallel to x with c 
pointing into the page. 


The integral over w is straightforward as it only in¬ 
volves exponential functions but it is very lengthy so it 
and the scalar product were performed symbolically using 
Matlab. Only the final inverse Fourier transform which 
gives the phase was evaluated numerically. 

A method to check the correctness of the solution is to 
plot the phase shift as a contour map for a = 90°. This 
gives the B-field projected through the thickness of the 
specimen and is shown in Fig. 01 It can be seen that the 
correct boundary conditions are fulfilled: the field lines 
spread as they approach the specimen surface from the 
interior and outside the specimen they are straight so 
that the field resembles that from a monopole if viewed 
far from the vortex. The figure shows that the effect of 
increasing the coherence length is to make the field 
less intense near the centre of the vortex as expected. 


D. Image Simulation 

Once the phase shift (j){x,y) has been calculated, the 
wavefunction of the electron beam is 'tpo{x,y) = 
and the intensity of the in-focus bright-held image is 
Io{x,y) = \tpQ{x,y)\'^. It is immediately clear that this 
gives I and an in-focus image is therefore featureless. In 
order to visualise flux vortices, out-of-focus images must 


where A is the electron wavelength. This is a convolu¬ 
tion and so is more conveniently evaluated as a multipli¬ 
cation in Fourier space via: 

= (44) 

After inverse transforming, the intensity in the out-of¬ 
focus image is given by /a/ = IV’A/P- 

III. EXPERIMENTAL METHOD 

MgB 2 single crystals were synthesised by the peritec- 
tic decomposition of MgNBg and their quality and bulk 
properties have been well characterised by a variety of ex¬ 
perimental techniquesi^^i^. The samples were prepared 
for electron microscopy at the Technical University of 
Denmark (DTU) using a Helios Nanolab focussed ion 
beam microscope (FIB). This is a dual-beam instrument 
in which a beam of gallium ions is used to mill the speci¬ 
men whilst secondary electrons emitted by the specimen 
are used to produce an image, however an electron beam 
can also be used to illuminate the specimen in order to 
take images without damaging the specimen. 

The MgB 2 single crystals were about I x I mm in the ab 
plane and about 100 |4.m thick in the c-direction. The in- 
situ lift-out technique was used to prepare the sample for 
electron microscopy. First, the FIB was used to deposit a 
3 |4m thick, 30 x 5 pm rectangle of platinum onto the ab 
surface to protect the sample beneath from ion damage. 
Trenches were then milled to a depth of 20 pm in the c 
direction around this to produce a slab standing in the 
centre of a crater. The top surface of the sample was 
smooth and this avoided the creation of the longitudinal 
thickness undulations reported in our last paper—. 



















A movable needle known as a micromanipulator was 
attached to the slab using platinum deposition and the 
slab was cut away from the rest of the specimen and 
extracted on the end of the micromanipulator. A sample 
can be tilted only to 25° in the electron microscope so 
to achieve a higher tilt angle, the FIB was used to cut a 
slot at 45° to the plane of an ‘Omniprobe’ grid. Further 
platinum deposition was used to attach the sample to 
this slot and the micromanipulator was then cut away, 
leaving the sample attached to the grid and tilted about 
its a axis by 45° with respect to the plane of the grid. 
The sample was then thinned to approximately 200 nm 
so that it was electron-transparent using a 30 kV Ga ion 
beam. Finally, the specimen surfaces were polished by a 
low-energy (2kV) Ga ion beam to minimise the damage 
layer caused by FIB milling. 

Electron microscopy was undertaken at DTU using an 
FEI Titan 80-300ST transmission electron microscope 
operated at 300 kV equipped with a Gatan imaging filter 
to record images. Under normal operating conditions, 
the main objective lens of the microscope applies a 2 T 
field to the specimen so to avoid this, the microscope was 
operated in low-magnification mode with the main objec¬ 
tive lens set to a low value and the image was focussed 
with the diffraction lens. Prior to imaging vortices, elec¬ 
tron diffraction was used to make a fine adjustment of a 
few degrees so that the tilt was purely about the a axis. 
Adjusting the sample so that it is tilted purely about a 
can be performed to better than 0.5° but this may alter 
the overall tilt angle and we judge that the tilt angle was 
a = 45 ±5°. 

The simulations were based on elastic electron scat¬ 
tering so experimental images were energy filtered so 
that only electrons which had lost O-IO eV on passing 
through the specimen contributed and an aperture was 
used so that only the 000 beam and the low-angle scat¬ 
tering from the vortices contributed to the image and 
the other crystallographic beams were excluded. The 
sample was cooled using a Gatan liquid-helium cooled 
‘IKHCHDTSOIO-Speciar tilt-rotate holder which has a 
base temperature of 10 K. 

The defocus and magnification were calibrated by ac¬ 
quiring images with the same lens settings as the original 
images from Agar Scientific’s ‘S106’ calibration specimen 
which consists of lines spaced by 463 nm ruled on an 
amorphous him. The defocus was found by taking dig¬ 
ital Eourier transforms of the images acquired from the 
calibration specimen and measuring the radii of the dark 
rings which result from the contrast transfer functioni^. 

A thickness map of the specimen was created by divid¬ 
ing an unhltered image by an energy-hltered image and 
taking the natural logarithm's which gives the thickness 
parallel to the electron beam, Z, as a multiple of the in¬ 
elastic mean free path, A^. To determine Xi, an electron 
hologram was taken at room temperature at an edge of 
the specimen which gives a phase shift proportional to 
the thickness, (p = CeVqI. Ce is a constant which de¬ 
pends only on the microscope voltage and has the value 


6.523 X 10® m~^V“^ at 300 kV. Vq, the mean inner po¬ 
tential, was calculated as Vq = 17.71 V from theoretical 
scattering factors given in ref. [U, giving = 244 ± 5 nm 
and the thickness, /, varied from 200-290 nm across the 
field of view of Eig. [5] Ideally the thickness of the whole 
specimen would have been determined by electron holog¬ 
raphy but the field of view of the interference region was 
not sufficiently large. 

A simplex algorithm^^ was used to minimise the re¬ 
duced value between the experimental images and 
simulations by fitting the vortex positions and the in¬ 
plane rotation angle of the vortices as well as Aab, Ac 
and The reduced value is defined as x^ = 

(1/A^) _ /simulation)2/^2 ^ jg ^j^g 

number of pixels used in the fit, Ij is the intensity of pixel 
j in the image and Uj is the noise associated with pixel 
j. We used (t| = jjaving previously taken a 

series of images of the vacuum with different electron in¬ 
tensities. A graph of the standard deviation versus the 
average intensity showed the noise was Shot noise (so 
that the square of the noise was proportional to the im¬ 
age intensity) and gave the proportionality constant c 
relating the counts recorded on the detector to the num¬ 
ber of electrons received. Unlike our previous publication 
where separate fits were made for each vortex^^, here all 
the vortex images were fit simultaneously with a single 
value of Aab, Ac and using the model described in 
section El 


IV. RESULTS 

Eig. [5Ka)-(f) show an experimental defocus series ac¬ 
quired at 10.8 K (the base temperature of our cooling 
stage) in a held of 4.8 mT. 68 images of vortices which 
were least affected by bend contrast were Ht and with 
a single value of A^b, Ac and along with the posi¬ 
tion of each vortex and its in-plane rotation angle using 
a simplex algorithm^^ to minimise the reduced x^ value. 
The specimen thickness and tilt angle were hxed at their 
calibrated values. 

Fig- Eig) shows the average of these images at each 
defocus level and (h) shows the average of the htted sim¬ 
ulations. To demonstrate that the ht is good, (i) shows 
the difference between images and simulations and (j) 
compares linescans taken across the vortex images. 

Aab, Ac and were then altered in turn and the er- 
rorbar on each judged by the point at which the differ¬ 
ence images displayed a discernibly worse ht as shown 
in Eig. Ela)-(d). This corresponded to an increase in 
the reduced x^ of 0.009 and the variation of x^ as each 
parameter is varies is shown in Fig. Ele). This yielded 
Aab = 107±33 nm, Ac = 134±6 nm and = 31±9 nm. 
The images are much more sensitive to the value of Ac 
than to Aab as a consequence of mounting the sam¬ 
ple tilted about the a axis: had it been tilted about 
c, the errors on Aab and Ac would be reversed. In a 
previous paper—, we obtained the more precise value of 
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FIG. 5. (a)-(f) Defocus series showing flux vortices in MgB 2 
at 10.8 K in a fleld of 4.8 mT with defoci: (a) A/ = —21.4, 
(b) —14.3, (c) —7.15, (d) 7.15, (e) 14.3 and (f) 21.4 mm. (g) 
Average vortex images at each defocus level, (h) Average sim¬ 
ulated images, (i) Difference images between experiment and 
simulation, (j) Linescans across the images from the region 
shown by the red box in (h). Black lines show experimental 
data and red lines the fit. 


Aab = 107 ± 8 nm using a sample cut in the ab plane. 

To check for amorphous ‘dead-layers’ of non¬ 
superconducting material on the sample surfaces caused 
by ion thinning, the thickness of the crystalline com¬ 
ponent of the sample was measured using the conver¬ 
gent beam diffraction technique described in ref. [H. 
This showed no difference between the total and crys¬ 
talline thicknesses to within the experimental error of 
±10 nm. Reducing the calibrated thickness values by 
50 nm to account for the largest conceivable dead-layer 
reduced Aab by 13 nm, reduced Ac by 27 nm and in¬ 
creased by 10 nm without altering the quality of the 
fit. Taking this into account gives: Aab = 100 ± 35 nm, 
Ac = 120 ± 15 nm and = 36± 10 nm. The anisotropy 
ratio in the penetration depth is then 7 a = 1.2 ± 0.4 and 
the coherence lengths are ^ab = = 41 ± 13 nm 

and ^c = = 34 ± 10 nm. Alternatively, us¬ 

ing the more precise value of Aab = 107 ± 8 nm and 
Ac = 120±15 nm gives 7 a = 1.12±0.16, ^ab = 39±11 nm 
and .^c = 35 ± 10 nm. 


V. DISCUSSION 

We now compare the values obtained here with those 
found using other techniques. These show that the con¬ 
ditions used (4.8 mT and 10.8 K) were in the low field 
limit (< 100 inT as established by neutron diffraction) 
but not quite in the low temperature limit (< 5 K from 
radio-frequency measurements). 

In 2005, Fletcher et al^ performed radio-frequency 
measurements which gave the change in the low-field pen¬ 
etration depth with temperature but not absolute values. 
At 10.8 K, Aab increased by 12±1 nm and Ac increased by 
18 ± 4 nm with respect to their low temperature values. 
Subtracting these from our most precise measurements 
of the penetration depths gives Aab = 95 ± 8 nm and 
Ac = 102 ± 15 nm and the anisotropy as 7 = 1.07 ± 0.18 
in the low field and low temperature limit. 

The most reliable measurement of the absolute value 
of the penetration depth is likely to be from neutron 
diffraction and in 2003 Cubitt et al^ found that at 
2 K, the extrapolated low field (< 100 mT) value was 
Aab = 82 ± 2 nm which is close to our value. 

As samples grow as thin plates in the ab plane, Cubitt 
et al. did not have direct access to Ac and so acquired 
diffraction patterns with vortices tilted at 45*^ with re¬ 
spect to the c-axis. It was uncertain whether the formula 
used to calculate the anisotropy was valid for a 2 -band 
superconductor— but data acquired at 2 K between 0.2- 
0.5 T indicated that 7 a varied with field and its extrap¬ 
olated value at low field was 7 a = 1-1 ± 0.3. In 2006, Pal 
et al^ used a different neutron diffraction technique to 
give 7 A = 1.1±0.1 at 4.9 K. Combining this with the 
neutron value for Aab gives Ac = 90 ± 8 nm which agrees 
with our value of Ac = 102 ± 15 nm. The anisotropy 
we obtain is close to the value of 1.01 calculated from 
first-principles by Golubov et al^ in the clean limit but 
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Penetration depth or coherence length/nm 


FIG. 6. Difference images between the average defocus se¬ 
ries and simulated images as each parameter is varied, (a) 
Difference images for the best-fit values of the parameters, 
(b) Difference images when Aab is varied: the top row shows 
difference images when Aab is set at two errorbars below the 
best-fit value. The next row is for Aab set one errorbar be¬ 
low. The next two rows are for Aab increased one and two 
errorbars above the best-fit value respectively, (c) A similar 
series showing the effect of changing Ac. (d) Series showing 
the effect of changing ^v- (e) Reduced values as Aab, Ac 
and are adjusted. Acceptable values of lie between the 
dashed lines. 


in common with Fletcher et al^, we find penetration 
depths approximately twice as large as predicted. 

Cubitt et al. interpreted their data assuming that the 
coherence length did not vary with field. If it did, the 
value they obtained, ^ab = 8± 1 nm, would apply only at 
high field (> 0.8 T). This is close to the value of 10 nm 
found from the upper critical field^^. 

Eskildsen et al^ measured the coherence length in 
2002 using scanning tunnelling microscopy (STM) to 
measure the width of vortex cores, scanning the ab plane 
with tunnelling in c. As the a carriers are confined to the 
a6-planes, the tunnelling current came almost exclusively 
from the tt band giving = 38.8 ± 0.7 nm at 0.32 K in 
a field of 50 mT (after adjusting for their slightly differ¬ 
ent model for the core). This agrees with our value of 

= 35 ± 10 nm at 10.8 K and 4.8 mT and supports this 
larger value at low held. 


VI. SUMMARY AND CONCLUSIONS 

We have described a new method to measure the pen¬ 
etration depth and coherence length of a superconductor 
in all directions at low applied magnetic held using trans¬ 
mission electron microscopy. The measurement does not 
need large-scale facilities and required one day to thin 
and mount the sample and another day to take the images 
required. The experiment was performed on a very small 
sample: 30 x 15 p.m and 200 nm thick so this method 
could prove useful for superconductors where only very 
small single crystals are available, as is the case for some 
iron-based superconductors. It is also useful if the pene¬ 
tration depth and coherence length vary with held, as is 
the case for MgB 2 , as the measurement is made at very 
low helds which can be difficult to access using other 
techniques. 

For a sample of MgB 2 cut in the ac plane and tilted to 
45 ±5° about the a-axis, we obtained Aab = 100 ±35 nm, 
Ac = 120±15 nm at 10.8 K in a held of 4.8 mT. The large 
error in Aab is a consequence of tilting the sample about 
the crystallographic a axis. Had it been tilted about c 
instead, the errors on Aab and Ac would be reversed. We 
obtained a more precise value of Aab = 107 ± 8 nm at 
10.8 K in our previous paper— in which the sample was cut 
in the ab plane. Using this value gives Aab = 107 ± 8 nm, 
Ac = 120 ± 15 nm, ^ab = 39± 11 nm and ^c = 35 ± 10 nm 
which agree well with measurements made using other 
techniques discussed in section IVl 

Obtaining the most precise values for the penetration 
depths and coherence lengths in all directions using this 
technique requires taking images of the vortex lattice 
with the sample tilted in more than one direction. It 
might be thought that the sample could be mounted in 
the plane of the support grid and tilted with the mi¬ 
croscope goniometer, hrst about a and then about c. 
However, the design of conventional liquid helium cooled 
holders does not allow tilting to an angle higher than 
a = 25°, so the sample was instead mounted to the sup- 












11 


port grid at 45 ± 5° to give sufficient contrast in the 
image. Thus, to investigate a new superconductor and 
obtain the most accurate measurement of the penetration 
depth and coherence length in all directions, it would be 
best to cut two samples and mount both to the grid, one 
tilted about a and the other tilted about c. 

We used the Ginzburg-Landau model for the magnetic 
structure of flux vortices but as MgB 2 is a two-band su¬ 
perconductor, the vortices may well have a more complex 
structure as described in ref. We obtained good fits 
to the vortex images and there was no indication of a 
more complex structure within the accuracy of the tech¬ 
nique but our simulation scheme would allow any model 


for the vortex structure to be used for image simulations 
provided that the magnetic vector potential for a vortex 
in an infinitely thick superconductor is known. 
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